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Abstract 

We extend a previous Landau-gauge study of subtractive renormalization of 
the fermion propagator Dyson-Schwinger equation (DSE) in strong-coupling, 
quenched QED4 to arbitrary covariant gauges. We use the fermion-photon 
proper vertex proposed by Curtis and Pennington with an additional correc- 
tion term included to compensate for the small gauge-dependence induced by 
the ultraviolet regulator. We discuss the chiral limit and the onset of dynam- 
ical chiral symmetry breaking in the presence of nonperturbative renormal- 
ization. We extract the critical coupling in several different gauges and find 
evidence of a small residual gauge-dependence in this quantity. 

I. INTRODUCTION 



Strong coupling QED 4 has been studied for some time within the Dyson-Schwinger equa- 
tion (DSE) formalism both for its intrinsic interest and also as the basis for abelianized 
models of nonperturbative phenomena in technicolor theories and QCD. For recent reviews 
of Dyson-Schwinger equations and their application and numerous references see for ex- 
ample Refs. The usual approach is to write the DSE for the fermion propagator or 
self-energy, possibly including equations for the photon vacuum polarization or the fermion- 
photon proper vertex. In a recent study M it was shown for the first time how to implement 
nonperturbative renormalization in a numerical way within the DSE formalism. In that 
work the calculations were carried out in quenched approximation in Landau gauge. Here 
we will extend these studies to arbitrary covariant gauges and will also discuss the chiral 
limit. 
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The DSE's are an infinite tower of coupled integral equations and so it is always necessary 
to truncate this tower at some point and introduce an Ansatz for any necessary undefined 
Green's functions. It is of course important to ensure that this Ansatz be consistent with 
all appropriate symmetries of the theory and that it have the correct perturbative limit. 
The resulting nonlinear integral equations are solved numerically in Euclidean space by 
iteration. Dynamical (or spontaneous) chiral symmetry breaking (DCSB) occurs when the 
fermion propagator develops a nonzero scalar self-energy in the absence of an explicit chiral 
symmetry breaking (ECSB) fermion mass. We will refer to coupling constants strong enough 
to induce DCSB as supercritical and those weaker are called subcritical. We write the 
fermion propagator as 

s(v) = z{p2) = 1 m 

W i>-M{j> 2 ) A{p 2 ) p 1 - B(p 2 ) y) 

where we refer to A(p 2 ) = 1/Z(p 2 ) as the finite momentum-dependent fermion renormal- 
ization and where M(p 2 ) = B(p 2 )/A(p 2 ) is the fermion mass function. In the massless 
theory (i.e., in the absence of an ECSB bare fermion mass) by definition DCSB occurs when 
M(p 2 ) ^ 0. 

Until relatively recently, most studies have used the bare vertex as an Ansatz for the one- 
particle irreducible (1-PI) vertex T u (k,p), H-O despite the fact that this violates the Ward- 
Takahashi Identity (WTI) [|Kj. The resulting fermion propagator is not gauge-covariant, 
i.e., physical quantities such as the critical coupling for dynamical symmetry breaking and 
the fermion mass pole are gauge-dependent pTl|— pL3|| . There have been several studies which 
attempted to make the fermion DSE gauge-covariant by using improved vertex forms which 
satisfy the WTI but which possess kinematic singularities in the limit of zero photon mo- 
mentum [l^JT^|. A general form for T u (k,p) which does satisfy the Ward Identity and which 
has no unphysical kinematic singularities was given by Ball and Chiu in 1980 ||15|| : it consists 



of a minimal longitudinally constrained term which satisfies the WTI, and a set of tensors 
spanning the subspace transverse to the photon momentum q. 

While the WTI is necessary for gauge-invariance, it is not a sufficient condition and in 
itself does not ensure gauge covariance of the fermion propagator. Furthermore, with many 
vertex Ansdtze the fermion propagator DSE is not multiplicatively renormalizable, which is 
equivalent to saying that overlapping logarithms are present. There has been much recent 
research on the use of the transverse parts of the vertex to ensure both gauge-covariant and 



multiplicatively renormalizable solutions |12],[T6|-g4[] , some of which will be discussed below. 



With the exception of Ref. [Q, studies have mostly neglected the issue of the subtrac- 
tive renormalization of the DSE for the fermion propagator. Typically these studies have 
assumed an initially massless theory and have renormalized at the ultraviolet cutoff of the 
loop integration, taking Z\ — Zi — 1. Where a nonzero bare mass has been used, it has 
simply been added to the scalar term in the propagator. While in some circumstances for 
the special case of Landau gauge this can be a reasonable approximation, it is in general 
incorrect. Although there have been earlier formal discussions of renormalization P, |T2| , |T9" 



the important step of subtractive renormalization had not been performed prior to the recent 
study in Landau gauge [|J. 

Here we present the results of a study of subtractive renormalization in the fermion DSE 
in arbitrary covariant gauge for quenched strong-coupling QED 4 . Note that here the term 
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"quenched" means that the bare photon propagator is used in the fermion self-energy DSE, 
so that Z 3 = 1 and there is no renormalization of the electron charge. This is a somewhat 
different usage to that found in lattice gauge theory studies, since in our study virtual 
fermion loops may still be present in the proper fermion-photon vertex. 

The organization of the paper is as follows: The formalism is discussed in Sec. [TT]. This 
section contains discussions of the DSE for the renormalized fermion propagator, the Ansdtze 
for the proper vertex, the subtractive renormalization procedure, the chiral limit, and renor- 
malization point transformations. Our detailed numerical results are presented in Sec. |U 
and we present our summary and conclusions in Sec. |TV|. 



II. FORMALISM 
A. Renormalized DSE 

The DSE for the renormalized fermion propagator, in an arbitrary covariant gauge, is 

S-\p) = Z 2 (ji, A)[/f-m (A)] -iZ 1 ( f i,A)e 2 ^ ^-^S{k)T v {k,p)D^{q) ; (2) 



(2nY 

here q = k — p is the photon momentum, /i is the renormalization point, and A is a regular- 
izing parameter (taken here to be an ultraviolet momentum cutoff). We write mo (A) for the 
regularization-parameter dependent bare mass. The renormalized charge is e (as opposed 
to the bare charge eo), and the general form for the renormalized photon propagator is 

with £ the renormalized covariant gauge parameter and £o = Z 3 (n,A)^ the corresponding 
bare one. Since we work in the quenched approximation, we have for the coupling strength 
and gauge parameter respectively a = e 2 /4ii = «o = eo/47r and £ = £0, and for the photon 
propagator we have 

D^(q) - D^(q) = j Lg»> + - ^1 1 . (4) 



B. Vertex Ansatz 



The requirement of gauge invariance in QED leads to a set of identities referred to as 
the Ward-Takahashi Identities (WTI). The WTI for the fermion-photon vertex is 

q^(k,p) = S- 1 (k)-S- 1 (p), (5) 

where q = k — p . This is a generalization of the original differential Ward identity, which 
expresses the effect of inserting a zero-momentum photon vertex into the fermion propagator, 
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dp u 



(6) 



The Ward identity Eq. (^j) follows immediately from the WTI of Eq. (|5p after setting to zero 
all but the v component of q, dividing both sides of the WTI by q u and then taking q u — ► 0. 
In general, for nonvanishing photon momentum q, only the longitudinal component of the 
proper vertex is constrained, i.e., the WTI provides no information on T^(k, p) = T^ u T v {p, k) 
for q ^ 0. [We use the notation 2>" = - {q fi q u /q 2 ) and CJ™ = (q^q u /q 2 ) for the 
transverse and longitudinal projectors respectively] In particular, the WTI guarantees the 
equality of the propagator and vertex renormalization constants, Z 2 = Z\ (at least in any 
reasonable subtraction scheme The WTI can be shown to be satisfied order- by-order 
in perturbation theory and can also be derived nonperturbatively. 

As discussed in fliyZoH, this can be thought of as just one of a set of six general requirements 
on the vertex: (i) the vertex must satisfy the WTI; (ii) it should contain no kinematic 
singularities; (iii) it should transform under charge conjugation (C), parity inversion (P), 
and time reversal (T) in the same way as the bare vertex, e.g., 



C- 1 T,(k,p)C = -Tl(-p,-k) 



(7) 



(where the superscript T indicates the transpose); (iv) it should reduce to the bare vertex in 
the weak-coupling limit; (v) it should ensure multiplicative renormalizability of the DSE in 
Eq. (H); (vi) the transverse part of the vertex should be specified to ensure gauge-covariance 
of the DSE. 



Ball and Chiu [15[ have given a description of the most general fermion-photon vertex 
that satisfies the WTI; it consists of a longitudinally- const rained (i.e., "Ball-Chiu") part 
rg C , which is a minimal solution of the WTI, and a basis set of eight transverse vectors 
T^(k,p), which span the hyperplane specified by C^T^ik.p) = (i.e., q u T"(k,p) = 0), 
where q = k — p. The minimal longitudinally constrained part of the vertex will be referred 
to as the Ball-Chiu vertex and is given by 



(k + pY 
k 2 — p 2 



[A{k 2 ] - A{p 2 )]- 



- \B{k 2 ) - B(p 2 )} 



(8) 



Note that since neither C^T^ik^) nor T^ v T^ c {k,p) vanish identically, the Ball-Chiu ver- 
tex has both longitudinal and transverse components. The transverse tensors can be conve- 
niently written as 



mk, P 

T?(k,p 
Tf(k,p 
T£(k,p 
If(k,p 

T?(k,p 
T£{k,p 



p ll (k ■ q) - k^{p ■ q) , 
\p»(k-q)-k^p-q)}(K 



^-qH: 



q 7 

-2 



fc"] - 2i(p 

+ (p + k)»4, 



kyk x p u 



0\v 



-iq^" , 
Y(P 2 ~ k 2 ) 

hp 2 - k 2 ) [y*(jrf + f) ~ P* 1 ~ k»] - % (k + pf k x p v a Xv 



lx k l 'p x a uX + k fl p'-p^^ 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 
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where we use the conventions g^ v = diag(l, — 1, — 1, — 1), { r )^-, r y l/ } = 2^", and a^ v = 
(i/2)[7 M , j v ]. Note that these tensors have been written in a different linear combination to 
the ones presented in Ref. M. A general vertex is then written as 

r»{k,p) = n c (k,p) + j2r l (k 2 ,p 2 ,q 2 )Tt(k,p) , (17) 

i=l 

where the r< are functions that must be chosen to give the correct C, P, and T invariance 
properties. 

Curtis and Pennington published a series of articles fT2||T9| -f2T| describing their specifi- 



cation of a particular transverse vertex term, in an attempt to produce gauge- covariant and 
multiplicatively renormalizable solutions to the DSE. In the framework of massless QED4, 
they eliminated the four transverse vectors which are Dirac-even and must generate a scalar 
term. By requiring that the vertex reduce to the leading log result for k ^> p they 

were led to eliminate all the transverse basis vectors except Tq, with a dynamic coefficient 
chosen to make the DSE multiplicatively renormalizable. This coefficient had the form 

r 6 (k 2 ,p 2 ,q 2 ) = - l -lA{k 2 )-A{p 2 )]/d{k,p) , (18) 

where d(k,p) is a symmetric, singularity-free function of k and p, with the limiting behavior 
lim fc 2 >p 2 d(k,p) = k 2 . [Here, A{p 2 ) = 1/Z{p 2 ) is their l/jFQo 2 ).] For purely massless QED, 
they found a suitable form, d(k,p) = {k 2 — p 2 ) 2 /{k 2 + p 2 ). This was generalized to the case 
with a dynamical mass M(p 2 ), to give 

d{Kp) = ^- P r + m + M^)? (19) 

They then showed that multiplicative renormalizability is retained up to next-to-leading- 
log order in the DCSB case. Subsequent papers established the form of the solutions for 
the renormalization and the mass |21| and studied the gauge-dependence of the solutions 



12| . Dong, Munczek and Roberts |22| subsequently showed that the lack of exact gauge- 
covariance of the solutions was due to the use of a momentum cutoff in the integral equations, 
since this type of regularization is not Poincare invariant. The fact that there is still some 
residual gauge- dependence in the physical observables such as the chiral critical point shows 
that with a momentum cutoff the C-P vertex Ansatz is not yet the ideal choice. Dong, 



Munczek and Roberts [E21 derived an Ansatz for the transverse vertex terms which satisfies 



the WTI and makes the fermion propagator gauge- covariant under hard momentum-cutoff 
regularization. 

Bashir and Pennington [23,^4| subsequently described two different vertex Ansdtze which 



make the fermion self-energy exactly gauge-covariant, in the sense that the critical point for 
the chiral phase transition is independent of gauge. Specific constraints they have assumed 



for the vertex are, in [23| that the transverse vertex parts vanish in the Landau gauge, and 
in that the anomalous dimension of the fermion mass function 7 m is exactly 1 at the 
critical coupling. Their work is a continuation of that of Dong, Munczek and Roberts, and 
indeed their vertex Ansatz corresponds to the general form suggested in p2[|. 

However, the kinematic factors T2 1 3,6,8 i n both vertex forms are rather complicated and 
depend upon a pair of as yet undetermined functions W\^{k 2 : p 2 ) which must be chosen to 
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guarantee that the weak-coupling limit of T M matches the perturbative result. Renormaliza- 
tion studies of the DSE using these new vertex Ansdtze should be interesting and represent 
a direction for further research. 

For the solutions to the fermion DSE using the C-P vertex, the critical point for the 
chiral phase transition has been shown to have a much weaker gauge-dependence than that 



for the DSE with the bare or minimal Ball-Chiu vertices 271. In this work we will use the 



Curtis-Pennington Ansatz as the basis for our calculations. 

The equations are separated into a Dirac-odd part describing the finite propagator renor- 
malization A(p 2 ), and a Dirac-even part for the scalar self-energy, by taking |Tr of the DSE 
multiplied by p'/p 2 and 1, respectively. The equations are solved in Euclidean space and so 
the volume integrals / d 4 k can be separated into angle integrals and an integral / dk 2 ; the 
angle integrals are easy to perform analytically, yielding the two equations which will be 
solved numerically. 

One refinement of our treatment of the C-P vertex in the present work is associated 
with subtleties in the ultraviolet regularization scheme. Although there have been some 



exploratory studies of dimensional regularization for the DSE ||28|| , this has not yet proven 
practical in nonperturbative field theory and momentum cutoffs for now remain the regular- 
ization scheme of choice in such studies. Naive imposition of a momentum cutoff destroys 
the gauge covariance of the DSE because the self-energy integral contains terms, related to 
the vertex WTI, which should vanish but which are nonzero when integrated under cutoff 



regularization [p2[p3f . In Appendix |A] we derive an expression for one such undesirable term 



and show how it may subtracted in a simple way from the regularized self-energy. We have 
also calculated some DSE solutions with the usual uncorrected UV cut-off method for com- 
parison purposes, but otherwise we use this "gauge-improved" regularization combined with 
the C-P vertex throughout this work. This will be commented on futher in the discussion 



of numerical results in Sec. Ill 



C. Subtractive Renormalization 

The subtractive renormalization of the fermion propagator DSE proceeds similarly to 
the one-loop renormalization of the propagator in QED. (This is discussed in and in ||29|| , 



p. 425ff.) One first determines a finite, regularized self-energy, which depends on both a 
regularization parameter and the renormalization point; then one performs a subtraction 
at the renormalization point, in order to define the renormalization parameters Zi, Z 2 , Z 3 
which give the full (renormalized) theory in terms of the regularized calculation. 

A review of the literature of DSEs in QED shows, however, that this step is never actually 
performed. Curtis and Pennington |12] for example, define their renormalization point at 
the UV cutoff. 

Many studies take Zi = Z 2 = 1 |T2| , p!6| , p!8| -|2l[] ; this is a reasonable approximation in 
Landau gauge in cases where the coupling a is sufficiently small (i.e., a ;$ 1), but if a is 
chosen large enough, the value of the dynamical mass at the renormalization point may be 



significantly large compared with its maximum in the infrared. For instance, in Ref. |12 
figures for the fermion mass are given with a = 0.97, 1.00, 1.15 and 2.00 in various gauges. 
For a = 2.00, the fermion mass at the cutoff is down by only an order of magnitude from 
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its limiting value in the infrared. In general for strong coupling and/or gauges other than 
Landau gauge this approximation is unreliable. 

As shown in Ref. [|J subtractive renormalization can be properly implemented in numer- 
ical DSE studies without such approximations. We begin with a summary of the renormal- 
ization procedure |1]|§]. One defines a regularized self-energy £'(/!, A;p), leading to the DSE 
for the renormalized fermion propagator, 

S-\p) = Z 2 (ji, A)[jf-mo(A)] - £'(/i, A;p) 

= i>- m(/i) - £(//; p) = A(p 2 ) $ - B{p 2 ) , (20) 

where the (regularized) self-energy is 

E'Gu, A; p) = iZ^, A)e 2 f t^T^O*! k)V\^- k, p)D Xu {^ (p - *)) ■ (21) 

[To avoid confusion we will follow Ref. |l| and in this section only we will denote regularized 
quantities with a prime and renormalized ones with a tilde, e.g. £'(/i, A; p) is the regularized 
self-energy depending on both the renormalization point fi and regularization parameter 
A and £(/i;p) is the renormalized self-energy] As suggested by the notation (i.e., the 
omission of the A-dependence) renormalized quantities must become independent of the 
regularization-parameter as the regularization is removed (i.e., as A — > oo). The self-energies 
are decomposed into Dirac and scalar parts, 

E'(/i, A; p) = £' d (/i, A; p 2 ) ff + £' s (/i, A; p 2 ) (22) 

(and similarly for the renormalized quantity, £(/z,p)). By imposing the renormalization 
boundary condition, 

S-\ P ) a _ a =y-m(Ai), (23) 

one gets the relations 

t d ^-p 2 ) = S^(/i, K-p 2 ) - E' d >, A;^ 2 ) (24) 

for the self-energy, 

Z 2 ( / x,A) = l + S' d (/i,A; /U 2 ) (25) 

for the renormalization constant, and 

m (A) = [m(ji) - S' s (/i, A; /i 2 )] /Z 2 (jjl, A) (26) 

for the bare mass. The mass renormalization constant is then given by 

Z m (ji,A) = mo{A)/m(ji), (27) 

i.e., as the ratio of the bare to renormalized mass. 

The vertex renormalization, Zi(p, A), is identical to Z 2 (/J,, A) as long as the vertex Ansatz 
satisfies the Ward Identity; this is how it is recovered for multiplication into £'(//, A; p) in 
Eq. (0). 
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In order to obtain numerical solutions, the final Minkowski-space integral equations are 
first rotated to Euclidean space. [Note that all equations in Sees. | and |T| are written in 
Minkowski space.] They are then solved by iteration on a logarithmic grid from an initial 
guess. The solutions are confirmed to be independent of the initial guess and are solved 
with a wide range of cutoffs (A), renormalization points (//), couplings (a), covariant gauge 
choices (£), and renormalized masses (m(/j,)). 

The chiral limit occurs by definition when the bare mass is taken to zero sufficiently 
rapidly as the regularization is removed. This is guaranteed, for example, by maintaining 
TOo(A) = as A — ► oo. Explicit chiral symmetry breaking (ECSB) occurs when the bare 
mass Too (A) is not zero, (or more precisely, whenever it is not taken to zero sufficiently 
rapidly as A — > oo). Dynamical chiral symmetry breaking (DCSB) is said to have occured 
when M(p 2 ) ^ in the absence of ECSB. As the coupling strength increases from zero, there 
is a transition to a DCSB phase at the critical coupling a c . Concisely, the absence of ECSB 
means that we can set toq(A) = 0, and the absence of both ECSB and DCSB (i.e., a < a c ) 
means that M(p 2 ), m(fi), and toq(A) simultaneously vanish. (Recall that in the notation 
that we use here, M(p 2 ) = B(p 2 )/A(p 2 ) and m(/x) = M(/i 2 ).) This is the same definition of 
the chiral limit that is used in nonperturbative studies of QCD, see e.g. Refs. [fl| |3|j30|] and 
references therein. Obviously, any limiting procedure where we take Too (A) — > sufficiently 
rapidly as A — ► oo will also lead to the chiral limit 0. 



D. Renormalization Point Transformations 



A renormalization point transformation is a change of renormalization scale [i.e., 1Z({1, //) 
for n — > //] such that the bare mass(es) and coupling(s) remain fixed for fixed regularization 
parameter (A) and fixed renormalization scheme. This ensures that the physical observables 
of the theory are invariant under such a transformation. This set of transformations is 
associative [7£(/x, , //') = 1Z(fjb, fj,")], contains the identity [7£(/z, /i) = X], and contains 

all inverses [1Z(fj,, /i') -1 = ^(/A/ 2 )] an d hence is called the renormalization group. 

For the purposes of the discussion here we will now explicitly indicate the choice of 
renormalization point by a /i-dependence of the renormalized quantities, i.e., A(fi;p 2 ) = 
1/Z(n;p 2 ), M(fi;p 2 ) = B (fi; p 2 ) / 'A{ji\ p 2 ) , etc. Note that Eq. @ implies that 



A(fi;p 2 ) 
B(fi-p 2 ) 



Z 2 (^,A)-^(/!,A;p 2 



E d (fx,A;p 2 



Z 2 (/i, A)toq(A) + S' s (/x, A;p 2 ) = m(n) + A;p 2 



(28) 



The renormalization point boundary condition in Eq. ( p3|) then leads to A; /z 2 ) = 0, or 
equivalently, to the two boundary conditions A(p,; fi 2 ) = 1 and M(/x; /z 2 ) = B(p,; fi 2 ) = to(^). 
From Eq. (p^) and the fact that Zi(/x, A) = Z 2 (fi, A) we have 



A(fi;p 2 )/Z 2 (ix,A) 
B(fi;p 2 )/Z 2 (^A) 



1 



V^A-p^/Z^A) 



m (A)+ E^A^/Z^A) 



(29) 



Consider the effects of an an arbitrary rescaling A(p 2 ) — > cA(p 2 ) and B(p 2 ) — ► cB(p 2 ), 
[i.e., M(p 2 ) fixed], for some real constant c. It is straightforward to see that under such 
a rescaling we have S(p) — > (l/c)5(p) and T"(p',p) — > cT u (p',p). It follows that the RH 
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sides of Eqs. (|29|) are unaffected by such an arbitrary rescaling. Hence, it follows that the 
choice of renormalization point boundary conditions is equivalent to the choice of scale for 
the functions A and B. 

Let us consider this observation in more detail. Since we are working in the quenched 
approximation, where e 2 and D are unaffected by a change of renormalization point, it follows 
from Eq. (|2l|) that A;p 2 )/Zi(/i, A) is renormalization point independent since a change 
of renormalization point is a rescaling of A and B. Then since mo (A) is renormalization 
point independent by definition, the entire RHS of Eqs. (|29"D must be independent of the 
choice of renormalization point. Thus, under a renormalization point transformation we 
must have for all p 2 

M( / i';p 2 ) = M( / i;p 2 ) = M(p 2 ), 
A^'-p 2 ) Z 2 (/z',A) , 1 

A(^p 2 ) Z 2 (»,A) A ^^ } A(fi;fi' 2 ) [M) 

from which it follows for the fermion propagator that S (/V '; p) / S (/i; p) = Z 2 (/i, A)/Z 2 (/i', A) 
in the usual way. The behavior in Eq. (|30|) is explicitly tested for our numerical solutions. 
It is clear from Eq. fl3~0|) that having a solution at one renormalization point (/j) completely 
determines the solution at any other renormalization point (//) without the need for any 
further computation. 

An alternative derivation of this result which starts from the renormalized action and 
which applies to the general unquenched case can be found for example in Sec. 2.1 of 
Ref. 0. For brevity we can denote the above renormalization point dependence of the 
fermion propagator by S(fjL]p) oc l/Z 2 (/x, A). In the general unquenched case pi we would 



have in addition D aL, (p,;q) cx oc l/Z 3 (/x, A), e(/i) cx Z 2 (/x, AWZ 3 (/x, A)/Zi(/x, A), and 



e(ji)V(ji- q,p) oc Z 2 (ji,A)JZ 3 ((jl,A). 



III. RESULTS 

Solutions were obtained in Euclidean space for the DSE for couplings a from 0.1 to 
1.30, in gauges with £ from -0.25 to 3, and with a variety of renormalization points and 
renormalized masses. All results in this section refer to Euclidean space quantities. In the 
graphs and tables that follow, there are no explicit mass units. Since the equations have 
no inherent mass-scale, the cutoff A, renormalization point /x, m(/z), and units of M(p 2 ) or 
B(p 2 ) all scale multiplicatively, and the units are arbitrary. In four dimensions the coupling 
has no mass dimension, therefore it remains unchanged for all such choices of mass units. 

Fig. ^ shows a family of solutions characterized by a = 1.00, fi 2 = 1 x 10 8 , m(/i) = 400, 
and gauge parameters from -0.25 to 1.25. We see that while A and B are strongly gauge 
dependent, the mass function M(p 2 ) = B(p 2 )/A(p 2 ) is relatively insensitive to £. The 
location of the mass pole of the physical electron must of course be independent of gauge, 
and this gauge independence has been demonstrated explicitly using the WTI for example by 



Atkinson and Fry [31]. Their proof assumes that the bare mass mo (A) is itself independent 
of gauge. Hence, in a fully gauge covariant treatment the mass function is independent of 
gauge at two scales (i.e., at the mass pole and at the UV regularization scale A). In our 
study we find that the mass function is relatively insensitive to the choice of gauge for all 
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p 2 . The nature of the Landau-Khalatnikov transformations makes the possibility of a 
gauge independent M(p 2 ) seem rather unlikely. 

The stability of the renormalized DSE solutions with respect to variations in the ultra- 
violet cutoff is evident in Fig. |2|. This graph shows solutions with a = 1.15, fi 2 = 10 8 , 
m(/i) = 400, and gauge £ = 0.25. The cutoff A 2 was varied over several orders of magni- 
tude with no apparent change in the solutions over the common range of momenta. This 
numerical stability was shown in other tests as well. For instance, we extracted the mass 
M{p 2 = 0) for solutions with a = 1.00, /i 2 = 10 4 , m(/x) = 0, and observed variations of less 
than one part in 10 4 as the UV cutoff was varied over six orders of magnitude. 

Tables |, |IT|, and [ITT] show the evolution of the renormalization constants Z 2 (fi,A), 
Z m (fi,A), and the cutoff-dependent bare mass m (A) as a function of the UV regulator 
A. We see that as we move further from Landau gauge ^(/x, A) decreases more rapidly with 
increasing A. In addition, we observe that the bare mass exhibits decaying oscillations with 
increasing A, which is directly related to the oscillations characteristic of the supercritical 
coupling and subtractive renormalization at large p 2 (see additional discussion later) . 

Fig. |3| shows solutions with the coupling varying from subcritical (a = 0.6) to supercrit- 
ical (a = 1.4) values, with identical renormalization point, renormalized mass and gauge. 
Here we see that the the nodes in the mass function M(p 2 ) move to lower momenta and the 
oscillations become more pronounced as the coupling is increased further and further above 
critical coupling. 

In order to test the gauge- invariance of the chiral critical point, we extracted the critical 
coupling from solutions in the Landau gauge and in two other covariant gauges, with £ = 0.25 
and 0.5. Miransky et ai, |||6| working in the quenched ladder approximation in Landau 
gauge, found that the infrared limit of the dynamical mass M(0) has an infinite-order phase 
transition, 



M(0) ~ 4A exp 



7T 



(31) 



Following their approach, we assume a similar form for the dynamical mass near the critical 
coupling, 



M(0) = Mexp 



(32) 



and construct an order parameter which is expected to have a second-order phase transition, 
— 1/ log[M(0)/M']. Since the inherent mass scale M is not known a priori, it is necessary to 
choose a reasonable scale M', and then calculate a corrected fit, which also yields the actual 
value of M. 

DSE solutions were obtained for several values of the coupling in each gauge, with the 
renormalization fi 2 = 10 4 . We set m(/z) = for the purpose of studying the transition. 
The IR mass limit was extrapolated for each solution and then the order parameters were 
calculated using an assumed mass scale M' = 200. The resulting critical curves are shown 
along with the M(0) values from the solutions, in Fig. |j. The parameters from the nonlinear 
fits are given in Table [TV]. The critical exponents j3 are the same to within their numerical 
tolerance, and suggest that (3 may be independent of gauge. Although the values of a c are 
close in value, there is clear evidence of residual gauge- dependence. 
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Our value for a c in the Landau gauge is very close to the value of 0.933667 found by 



Atkinson et al |27j in their bifurcation analysis of the solutions of the fermion DSE with 
the Curtis-Pennington vertex. In addition we find that the critical coupling varies with £ 
in the same direction. Recall that these authors used the unrenormalized equations and 
relaxed the UV momentum cutoff to infinity in order to remove cutoff artifacts, whereas we 
have used subtractive renormalization and our gauge covariance correction. Hence we can 
anticipate a small difference between the critical couplings in our approaches. 

Fig. |5| shows a family of equivalent solutions renormalized at different momentum scales 
fi and these provide a direct check on the behavior predicted in Eq. (PU|). All have coupling 
a = 1.15 and £ = 0.5, and the renormalization scale fi 2 is stepped by powers of 10, with 
renormalized masses m(/z) chosen such that m(/i) = M (/x 2 ) for each //. It is clear that the 
resulting mass curves are identical at all p 2 as expected, and in each case A(/j,;p 2 ) scales as 
predicted in Eq. ( |3~0"D showing that the renormalized DSE does transform correctly. 

The anomalous dimension of the mass, j m , is defined by the asymptotic scaling of the 
dynamical mass with p 2 , 

9\ (7m/2)"l 



M(p 2 ) ~ (33) 

As well as depending on the coupling, it shows a slight dependence on the gauge, as shown 
in Fig. |6|. The dynamical masses shown are from DSE solutions with a = 0.5, fi 2 = 10 8 , 
and m(/z) = 400, in Landau gauge and in gauges £ = 0.25 and 0.5. They are scaled by 
multiplication with (p 2 / 'yU 2 ) 1 ^ 7 " 1 / 2 ), where the value of y m used was that extracted from 
the Landau gauge solution, y m = 1.716638. The gauge-dependence of y m shows up as the 
slight difference in slopes on the log- log plot. (The dips apparent at the end of the curves 
are due to having a hard momentum cut-off A 2 = 10 16 . As A is increased these move to 
higher momenta also.) In Landau gauge, the actual power of 1/p 2 with which M(p 2 ) falls 
asymptotically is s = l — (y m /2) = 0.141681. For the gauges £ = 0.25 and 0.5, the anomalous 
dimensions are 1.713948 and 1.711274, giving powers of 1/p 2 equal to 0.143026 and 0.144363 
respectively. 

Miransky has studied the form of the mass renormalization Z m in the bare vertex approx- 
imation in Landau gauge, and without subtractive renormalization |7||§. In this treatment 
he finds Z m {ji, A) = (/x 2 /A 2 ) 2~ 7 ', with the exponent 



1 '(a) = -^/l-a/a c , (34) 

where the critical coupling for DCSB in that approximation is a c = ir/3. This would imply 
an asymptotic scaling for the dynamical mass that goes like M(p 2 ) ~ (p 2 ) 7 ~ 1//2 , so that 
the anomalous mass dimension would be related to 7' by 7' = (y m — l)/2. Recent articles 
by Holdom [^] and Mahanta [51,S5| claim that for quenched theories at criticality the 



mass anomalous dimension y m should be exactly 1, giving M(p 2 ) ~ 1/p as in the bare 
approximation, and that in particular this result should be independent of the gauge ||35|| . 



In Miransky's treatment this corresponds to the vanishing of 7' at the critical coupling. We 
find for subcritical couplings that 7 m ~ 1.71 and that it is gauge insensitive. Hence, as for 
the pure CP- vertex case, our gauge corrected CP vertex does not lead to y m = 1. 
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The gauge covariance correction described in the Appendix leads to an exact restoration 
of gauge covariance in the subcritical case with no explicit chiral symmetry breaking, i.e., 
no bare mass (see, e.g., Ref. [^2j). The difference between a DSE solution with naive cutoff 
regularization and those with the gauge covariance correction is shown in Fig. [7]. Both 
solutions have a = 1.15, £ = 0.5, and are renormalized at fi 2 = 10 s and m(/i) = 400. The 
quantitative change induced by our gauge covariance correction was found to be relatively 
small in the presence of a substantial mass function M(p 2 ). 

We find, as in our previous study in Landau gauge ||, that for supercritical couplings 
the dynamic mass crosses zero; for solutions away from the chiral limit, the position of the 
first node depends on the gauge, as shown in Fig. In fact as the cutoff is increased M(p 2 ) 
shows damped oscillations periodic in logp 2 , as shown in Fig. ||. This has been discussed by 
several authors |7| J35|J36|1 ; in particular it is shown using some simplifying approximations, 



m 



3q] that in Landau gauge, the DSE reduces to a differential equation for M{p) which has 



the solution 

M{p) = k (£\ cos Qln(p 2 //i 2 ) v / «K- 1 + ^) > (35) 

with n cos cj) = m(fi). However, the approximations used in deriving this result are not 
applicable outside Landau gauge and even in Landau gauge lead to differences from the 
present treatment. We find that the functional form is substantially correct, but that both 
the mass dimension and the period of oscillations depend on the coupling and the choice 
of gauge. In fact, for the case shown in Fig. || the mass dimension that fits |M(p 2 )| is 
7 m = 1.115. The dependence of the period on a and a c is also not as simple as that in Eq. 



IV. SUMMARY AND CONCLUSIONS 

We have extended our previous work on the numerical renormalization of the DSE |4j] to 
arbitrary covariant gauges. The procedure is straightforward to implement and extremely 
stable. It becomes numerically more challenging for covariant gauges far removed from 
Landau gauge and for large couplings (a 1). The importance of the approach is that it 
removes the issue of cut-off dependence and allows solutions to be obtained for any choice 
of renormalization point. We have described the procedure for performing renormalization 
group transformations between solutions with different renormalization points. We saw that 
a knowledge of the solution at one renormalization point automatically provides the solutions 
at all renormalization points. 

This then allows also comparisons with results from lattice studies of QED, which should 
prove useful in providing further guidance in the choice of reasonable Ansatze for the vertex 
and photon propagator. Without renormalization only the unrenormalized, regulated quan- 
tities would be be obtained and any such comparisons would be meaningless. In addition, in 
order to study the nonperturbative behaviour of renormalization constants such as Z\(fi, A), 
Z 2 (/U,A), and Z m (//, A) they must be numerically extracted and so a method such at that 
described here would be essential. 

The context of this study has been quenched four-dimensional QED with a modified 
Curtis-Pennington vertex, since that vertex Ansatz has the desirable properties of making 
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the solutions approximately gauge-invariant and also multiplicatively renormalizable up to 
next-to-leading log order. The technique described can be generalized to apply elsewhere 
(e.g., QCD), whenever numerical renormalization is required. 

The solutions are stable and the renormalized quantities become independent of regular- 
ization as the regularization is removed, which is as expected. For example, the mass function 
M(p 2 ) and the momentum-dependent renormalization A(p 2 ) = 1/Z(p 2 ) are unchanged to 
within the numerical accuracy of the computation as the integration cutoff is increased by 
many orders of magnitude. The mass renormalization constant Z m (ji,A) converges to zero 
with increasing A because the mass function M(p 2 ) falls to zero sufficiently rapidly at large 
p 2 . The absence of divergences of Zi(fi, A) = Z 2 (/x, A), m (A), and Z m (/i,A) in the limit 
A — > oo is a purely nonperturbative result and is in sharp contrast to the perturbative case 
where these constants diverge at all orders. 

In order to study the critical point and exponents for the transition to DCSB, one sets 
the renormalized mass m(/i) to and varies the coupling. For subcritical couplings m (A) 
remains zero, and the dynamical mass M(p 2 ) is identically zero, while for supercritical 
couplings M(p 2 ) ^ [and mo (A) ^ for finite A]. We have extracted the critical coupling 
for DCSB in Landau gauge (£ = 0), and in gauges with £ = 0.25 and 0.5. Our Landau gauge 
result is very close to the value found in P71; the values in the other gauges show a small 



residual gauge-dependence. 

For subcritical couplings, we find that the mass renormalization Z m (/x, A) scales approx- 
imately as Z m (ji,A) oc (ii 2 /A 2 f-^ a ^/ 2 , where e.g. 7 m (0.5,0) = 1.716638. Above the crit- 
ical coupling, the mass function shows damped oscillations around zero, periodic in ln(p 2 ). 
We have extracted mass anomalous dimensions 7 m for some subcritical and supercritical 
couplings, and find them all greater than 1. 

We have shown that our modified Curtis-Pennington vertex, while removing the viola- 
tion of gauge covariance in the massless subcritical case (i.e., when there is no explicit or 
dynamical chiral symmetry breaking), has not been sufficient to remove the small residual 
violation of gauge covariance in the general case. Hence it is important to attempt to extend 
this work to include other regularization schemes (e.g., dimensional regularization) and to 
vertices of the Bashir- Pennington type [0,0 ■ 
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APPENDIX A: UV REGULATOR AND GAUGE COVARIANCE 



Regulators applied to divergent integrals in field theory always destroy some continuous 
symmetry, and in particular the use of a momentum cutoff destroys gauge covariance. This 
appendix describes a modification of the self-energy integrals in the regularized DSE, which 
will at least partially restore this symmetry. 

The basis of this change in the regularization scheme is that when the self-energy E, given 
in Eq. fl2lD , is evaluated under cutoff regularization, it contains a term related to the vertex 
WTI, which should vanish but which integrates to give a nonzero contribution because the 



cutoff regularization scheme is not translationally invariant p2| , p3[ . We therefore evaluate 
this term separately and subtract it from the self-energy. It turns out to affect only the 
value of the renormalization Z 2 (/i, A). 

We write the photon propagator, in Minkowski momentum space and in an arbitrary 
covariant gauge, as 



-sv + 



1 



1 1 



q 2 J 1 + II(g 2 ) q 2 s q 2 q 2 

dim + dim- 



(Al) 

(A2) 
(A3) 



where T^viq) and C^q) are the transverse and longitudinal projectors respectively, T^ v (q) = 

g»v - q^/q 2 , £»M = Wvlq 2 - 

The renormalized fermion propagator is as in Eq. (|l|). T u (k,p) is the renormalized proper 
vertex; q = k — p is the photon momentum. The renormalization point for S, T, and D is 
p 2 = fi 2 , however, we will not always write it explicitly. 

The "naively" regularized self-energy (under a regularization scheme with parameter A) 



is 



S A (p) = iZ 1 ( f ,, A)e 2 f ^L^S(k)T»(k,p)D,M 



(2tt)< 



(A4) 



(where we write the regularization as though it were a momentum cutoff, but it need not 
be). If T u (k,p) satisfies the WTI, q u T u {k,p) = S~ 1 (k) — S^ 1 (p) , we can rewrite the DSE as 



S A (/i,A;p) = zZ 1 (/i,A)e 2 y' 



A ^YS(k)T»(k,p)Dl(q) 



(2tt)< 



q 2 q 2 



(A5) 



^(/i,A)e 2 f ^rS(k)T"(k,p)DlM 



iCZx(ii,A)e 



., , A d 4 k i\ 



(2vr) 4 q 2 q 2 



+iiZ x {^ A)e 2 J 



(A6) 
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The boxed integral is odd in q, and should vanish in any translationally invariant regu- 
larization scheme; otherwise, it contributes and we expect it to destroy the gauge-covariance 
of S A . However, we can define a "gauge- improved" self-energy by cancelling this undesirable 
term 

£'(//, A; p) = £ V A; P) + A)e 2 f -0- ^1 (A7) 

and since the added integral is Dirac-odd, upon decomposing into scalar and spinor parts 
as in Eq. (|22| ) we have 

S;( / i,A;p 2 ) = S^( yU ,A;p 2 ), (A8) 

£',(//, A; p 2 ) = S^(/i,A;p 2 ) +<eZ 1 ( A i, A)e 2 / A ^^| . (A9) 

This modification combined with the C-P vertex is the Ansatz used in all calculations in 
this work unless explicitly stated otherwise. 

Converting to Euclidean metric (but suppressing the "Euclidean" subscript on momenta 
for convenience) gives 

£^,A;p 2 ) = S^(/i,A;p 2 ) -^i(^A)e 2 / A ^^| . (A10) 

Since in our case the regularization is an ultraviolet cutoff in momentum, we have 

£',(/,, A; p 2 ) = EjfoAjp 8 ) -»,A)^ k 2 dk 2 J\m 2 9d9^ . (All) 



Introducing the variables x = p 2 , y = k 2 , z = (k—p) 2 = x + y — 2^/xy cos9, x > = m&x(x,y), 
x< = min(x,y), we have 



xz 2 



^,A;P 2 ) = ^,A;p 2 )-^,A)^- 2 f A ydy [\m 2 9d0^ C ° s9 X 

2ix z J Jo 

= E d (/i,A;p ) -Z x (/i,A)— - - Jf 



27T 2 2 x Jo (x> — x<)x> 
E^A^ + Z^A)^, (A12) 



57T 



Thus the gauge covariance correction that we use here is to cancel the boxed term in Eq. ( |A6| ) 
by adding Z 1 (fi, A)a£/87r to the naive regularized self-energy in Eq. (|A4j) and this is how 
it is implemented in our program. The effects of including and not including this gauge 
covariance correction can be studied numerically. 
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FIG. 1. The finite renormalization A(p 2 ) and the mass function M(p 2 ) are shown for various 
gauge parameters £. These results have coupling a = 1.00, renormalization point /i 2 = 10 8 , and 
renormalized mass m(fj,) = 400. In the low p 2 region the larger gauge parameter has the larger 
value of M(p 2 ). 
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FIG. 2. The finite renormalization A(p 2 ) and the mass function M(p 2 ) are shown for various 
choices of the regularization parameter (i.e., ultraviolet cut-off) A. These results have coupling 
a = 1.15, renormalization point /x 2 = 10 8 , renormalized mass m(/i) = 400, and gauge parameter 
£ = 0.25. The stability of the subtractive renormalization procedure is apparent. 
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FIG. 3. The finite renormalization A(p 2 ) and the mass function M{jp 2 ) are shown for various 
choices of the coupling strength a. These results have renormalization point fj, 2 = 10 8 , renormalized 
mass m(n) = 400, and gauge parameter £ = 0.50. 
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Critical Curves for gauges t, = 0, 0.25, 0.5 




0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 

Coupling a 

FIG. 4. The critical curves for three choices of gauge parameter showing the existence of 
residual gauge-dependence in the Curtis-Pennington vertex. For the purposes of studying the chiral 
phase transition all solutions were renormalized with the choice, m(fi) = 0, with the renormalization 
point /i 2 = 10 4 . The order parameter is evaluated using an arbitrary reference mass scale choice of 
M' = 200.0. Diamonds (O) connected by the solid smooth curve, are order parameter values for 
the Landau gauge; pluses (+) connected by the dashed smooth curve, are values for £ = 0.25; and 
boxes (□) connected by the dot-dashed smooth curve, are values for £ = 0.5. 
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FIG. 5. The finite renormalization A(p 2 ) and the mass function M{p 2 ) are shown for various 
choices of renormalization point. These results have coupling strength a = 1.15 and gauge param- 



eter £ = 0.50. Each of these results corresponds to M(p 2 ) = 400 at p 2 = 
renormalization point independent and A(p 2 ) varies as described in Eq. pC 



10 8 . Hence, M{p 2 ) is 
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FIG. 6. Asymptotic mass scaling below critical coupling as a function of gauge, with a = 0.5, 
fj, 2 = 10 s , m(fi) = 400. The scaling applied to the masses uses the anomalous dimension found for 
the Landau guage, 7 m (£ = 0) = 1.716638. the extracted anomalous dimensions for the other two 
curves are 7 m (0.25) = 1.713498, and 7 m (0.5) = 1.711274. 
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FIG. 7. The finite renormalization A(p 2 ) and the mass function M(p 2 ) are shown with and 
without the gauge covariance correction. We see that the correction is a relatively small effect. 
These results have coupling strength a = 1.15, renormalization point fi 2 = 10 8 , renormalized mass 
m{n) = 400, and gauge parameter £ = 0.50. 
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Detail of the node in the mass function M(p 2 ) for various gauge choices. These 
coupling strength a = 1.15, renormalization point /J, 2 = 10 8 , and renormalized mass 
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FIG. 9. Absolute value of the dynamical mass, showing damped oscillations periodic in ln(p 2 ). 
The solution shown has a = 1.25, £ = 0.25, and is renormalized with pL 2 = 10 4 , m(/x) = 0. 
The powerdaw fit, which runs tangent to the dynamical mass curve, is C(p 2 / y 2 )^™/ 2 )- 1 , with 
C = 4.394 x 10~ 2 , 7m = 1.115. 
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TABLES 



TABLE I. Renormalization constant ^(/i, A), bare masses mo (A), and mass renormalization 
Z m (fi,A), as a function of UV cutoff for a = 1.15 in the Landau gauge (£ = 0). All solutions are 
with renormalization point fj? = 1.00 x 10 8 and renormalized mass m(/x) = 400.0 



A 2 


^ 2 (/i,A) 


m (A) 


Z m (/j,,A) 


1 x 10 8 


0.9999135 


2.306 x 10 2 


5.765 x lO" 1 


1 x 10 9 


0.9998483 


5.358 x 10 1 


1.339 x lO" 1 


1 x 10 10 


0.9998468 


4.443 


1.111 x 10~ 2 


1 x 10 11 


0.9998469 


-3.932 


-9.831 x 10~ 3 


1 x 10 12 


0.9998469 


-2.847 


-7.117 x 10~ 3 


1 x 10 13 


0.9998469 


-1.182 


-2.954 x 10~ 3 


1 x 10 14 


0.9998469 


-3.408 x lO" 1 


-8.520 x 10~ 4 


1 x 10 15 


0.9998469 


-5.390 x 10~ 2 


-1.348 x 10~ 4 


1 x 10 16 


0.9998469 


1.043 x 10~ 2 


2.607 x 10~ 5 


1 x 10 17 


0.9998469 


1.276 x 10~ 2 


3.191 x 10~ 5 


1 x 10 18 


0.9998469 


6.171 x 10~ 3 


1.543 x 10~ 5 


1 x 10 19 


0.9998469 


2.042 x 10~ 3 


5.105 x 10~ 6 



TABLE II. 


Renormalization constant ^(/i, A), bare masses mo (A), and mass 


renormalization 


Z m (n,A), as a 


function of UV cutoff for a 


= 1.15 in the gauge with £ = 0.25. 


All solutions are 


with renormalization point /i 2 = 1.00 x 10 8 


and renormalized mass m(/z) = 400.0 




A 2 


Z 2 (V,A) 


m (A) 


Z m ([i,A) 


1 x 10 8 


0.999943 


2.239 x 10 2 


5.598 x 10" 1 


1 x 10 9 


0.9486 


4.918 x 10 1 


1.229 x 10" 1 


1 x 10 10 


0.8999 


2.034 


5.085 x 10~ 3 


1 x 10 11 


0.8537 


-4.898 


-1.225 x 10~ 2 


1 x 10 12 


0.8099 


-3.102 


-7.756 x 10- 3 


1 x 10 13 


0.7683 


-1.193 


-2.981 x 10 -3 


1 x 10 14 


0.7289 


-3.059 x 10- 1 


-7.647 x 10- 4 


1 x 10 15 


0.6915 


-2.886 x 10~ 2 


-7.214 x 10~ 5 


1 x 10 16 


0.6560 


2.145 x 10~ 2 


5.362 x 10~ 5 


1 x 10 17 


0.6224 


1.609 x 10~ 2 


4.023 x 10" 5 


1 x 10 18 


0.5904 


6.679 x 10~ 3 


1.670 x 10~ 5 
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TABLE III. Renormalization constant Ziiji, A), bare masses mo (A), and mass renormalization 
Z m {fjL, A), as a function of UV cutoff for a = 1.15 in the gauge with £ = 0.5. All solutions are with 
renormalization point \j? = 1.00 x 10 8 and renormalized mass m(/x) = 400.0 



A 2 


^ 2 (M,A) 


mo (A) 


Z m (p,,A) 


1 x 10 8 


0.99997 


2.176 x 10 2 


5.441 x 10" 1 


1 x 10 9 


0.8999 


4.513 x 10 1 


1.128 x 10 _1 


1 x 10 10 


0.8099 


-1.455 x lO" 1 


-3.638 x 10~ 4 


1 x 10 11 


0.7289 


-5.736 


-1.434 x 10~ 2 


1 x 10 12 


0.6560 


-3.299 


-8.248 x 10~ 3 


1 x 10 13 


0.5904 


-1.181 


-2.954 x 10~ 3 


1 x 10 14 


0.5314 


-2.653 x lO^ 1 


-6.632 x 10~ 4 


1 x 10 15 


0.4783 


-3.676 x 10~ 3 


-9.190 x 10~ 6 


1 x 10 16 


0.4304 


3.151 x 10~ 2 


7.877 x 10~ 5 


1 x 10 17 


0.3874 


1.870 x 10~ 2 


4.674 x 10~ 5 



TABLE IV. Critical parameters for three choices of gauge, £ = 0, 0.25, and 0.50. These are 
extracted from nonlinear fits to the data in Fig. ||, using the form in Eq. (p2|). 



Parameter 


Landau (£ = 0) 


£ = 0.25 


£ = 0.5 


c 


2.877 ± .027 


2.858 ± .043 


2.851 ± .055 


a c 


0.93307 ± .00023 


0.92076 ± .00048 


0.90946 ± .00071 


P 


0.512 ± .003 


0.514 ± .005 


0.516 ± .007 


M 


154.3 ±5.2 


148.5 ±7.7 


145.4 ±9.4 


x 2 /dof 


0.0959 


.0388 


.0211 
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